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ABSTRACT 

Natural  circulation  loop  works  on  the  basic  principle  of  thermosiphon  effect,  which  doesn ’t  require  any  external  force  to 
drive  the  fluid  in  a  loop.  Hence  they  are  highly  applicable  in  power  industries  where  there  are  a  source  and  a  sink. 
However,  the  flow  pattern  in  a  single-phase  natural  circulation  loop  is  quite  unstable  especially  the  horizontal  heater 
horizontal  cooler  configuration  which  results  in  lower  loop  effectiveness.  A  lot  of  research  has  been  carried  out  to  study 
these  instabilities  and  find  a  way  to  curb  them.  Techniques  such  as  tilting  of  loop,  use  of  nanofluids,  adding  flow 
resistances  etc.  have  shown  promising  results.  But  they  have  their  own  practical  constraints.  Hence  in  the  present 
experimental  study,  a  Tesla-type  passive  valve  is  incorporated  in  natural  circulation  loop.  Due  to  its  presence,  not  only 
stability  period  is  reduced  (at  least  90  s  compared  to  non-Tesla  loop )  but  also  unidirectional  flow  is  achieved.  The  Tesla 
valve  performed  better  for  a  broad  range  of  heater  power  (90-170W)  and  was  also  able  to  provide  noticeable  resistance  to 
the  flow  when  the  loop  was  inclined  by  2.5°  in  the  opposite  direction. 
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1.  INTRODUCTION 

The  term  circulation  refers  to  a  closed-circuit  carrying  fluid  that  will  be  under  motion.  Circulation  or  flow  can  be  of 
two  types  viz.  natural  and  forced.  The  forced  type  is  catered  by  pumps,  blowers  etc.  based  on  the  application  and 
fluid  to  be  transported.  The  natural  circulation  (NC)  also  referred  to  as  thermosiphon  or  natural  convection  systems 
which  do  not  require  any  physical  component  to  trigger  its  motion  but  needs  a  channel  or  path  for  it  to  travel.  The 
driving  force  for  the  flow  occurs  due  to  the  presence  of  heat  source  and  heat  sink.  This  is  the  reason  why  natural 
circulation  loop  (NCL)  finds  application  in  power  industries,  solar  water  heaters,  transformer  cooling,  geothermal 
power  extraction,  cooling  of  internal  combustion  engines,  gas  turbine  blades  and  nuclear  reactor  cores. 

In  heat  pipe,  the  heating  and  cooling  happen  in  the  same  leg/channel  whereas,  in  NCL  there  are  two  separate 
legs  in  which  fluid  can  take  either  direction  once  getting  heated.  NCL  behaves  differently  for  the  various  heater  and  cooler 
orientations  as  thoroughly  researched  by  both  numerically  and  experimentally.  Horizontal  heater  horizontal  cooler 
(HHHC)  is  the  least  stable  arrangement  (Figure  1)  and  has  been  researched  a  lot  in  order  to  improve  its  directional 
stability  and  reduce  the  time  it  takes  to  form  stable  flow.  The  fluid  in  the  horizontal  position  can  take  either  of  the 
clockwise  or  anti-clockwise  direction,  giving  rise  to  uncertainty  and  ineffectiveness.  However,  vertical  heater  vertical 
cooler  (VHVC)  is  the  most  stable  configuration  as  flow  is  readily  established  forming  a  unidirectional  and  stable  loop. 
Other  orientations  are  horizontal  heater  vertical  cooler  (HHVC)  and  vertical  heater  horizontal  cooler  (VHHC). 
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Due  to  the  lack  of  control  and  interplay  of  various  forces  in  NCL,  various  kinds  of  instability  can  arise.  The  flow 
is  considered  stable  when  the  parameters  like  velocity  and  temperature  do  not  change  with  time.  All  NC  systems  take  a 
certain  amount  of  time  to  stabilize.  The  alternating  effects  of  buoyancy  and  gravity  due  to  change  in  density  across  the  loop 
causes  such  oscillation.  It  also  includes  the  actions  of  friction  and  velocity.  The  flow  in  HHHC  type  loop  may  take  any 
direction,  and  to  have  a  more  stable  and  controllable  system,  flow  need  to  be  unidirectional. 

2.  TESLA  VALVE 

A  Tesla  valve  is  a  valvular  conduit,  fixed-geometry  passive  check  valve  without  moving  parts.  It  allows  fluid  to  flow 
preferentially  in  one  direction.  The  device  is  named  after  Nikola  Tesla,  who  was  awarded  a  patent  in  1920  for  its  invention 
[2].  To  increase  the  efficiency  of  the  Tesla  valve,  one  has  to  increase  the  resistance  in  the  reverse  side  and  also  have  a 
lower  pressure  drop  in  the  forward  direction.  In  an  attempt  to  obtain  this,  D-type  Tesla  valve  is  designed  as  shown  in 
Figure  2.  However  the  old  version  is  depicted  in  Figure  3. 

3.  LITERATURE  REVIEW 


In  this  section,  the  recent  literature  review  related  to  single-phase  HHHC  type  of  NCL,  instabilities  and  restraining 
techniques  in  NCL,  and  Tesla  valve  have  been  discussed. 


To  understand  the  fluid  behaviour,  Misale  et  al.  [4]  have  carried  out  experiments  in  a  rectangular  loop  filled  with 
water  and  water-glycerol  mixture.  It  was  found  that  the  temperature  variations  throughout  the  loop  followed  a  pattern 
of  damped  oscillations  until  it  reached  a  steady  state.  The  one-dimensional  analytical  procedure  by  Keller  [5]  for  a  loop 
indicated  that  the  fluid  took  a  particular  direction  and  exhibited  periodic  motion  in  terms  of  its  magnitude  of  velocity. 
These  oscillations  were  observed  to  be  mainly  dependent  on  the  frictional  and  buoyancy  forces.  Vijayan  [6]  has 
experimentally  captured  the  characteristics  of  steady-state  and  stability  behaviour  of  NCL.  It  was  verified  that  the  steady- 
state  Re  in  uniform  or  non-uniform  diameter  NCL  could  be  expressed  as, 


Figure  1:  NCL-  HHHC  System  [1]. 


Figure  3:  Cross-Section  of  a  Tesla  Valve  [2]. 
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Ress=  C[  Gr„/NGf  (1) 

Where  the  constants  r  and  C  depend  on  the  nature  of  the  flow  (laminar  or  turbulent),  NG  is  the  geometric  ratio  and 
Grm  is  the  modified  Grashoff  number.  Further,  the  stability  criteria  for  uniform  and  non-uniform  diameter  loop  was  defined  as 
a  function  of  Stanton  number,  modified  Grashoff  number  and  geometric  ratio.  The  effect  of  wall  friction  on  stability  of  NCL 
in  different  flow  regimes  was  numerically  studied  by  Ambrosini  et  al.  [7]  and  concluded  that  the  friction  model  to  be  used  in  a 
loop  is  very  much  dependent  on  the  geometry.  Swapnalee  and  Vijayan  [8]  have  done  analytical  and  experimental  studies  with 
respect  to  four  different  orientations  of  the  heater  and  cooler  and,  presented  general  flow  equations  assuming  partial  laminar 
and  partial  turbulent  flow.  The  obtained  equations  were  in  good  agreement  with  the  data  in  earlier  literary  works. 

Yadav  et  al.  [9]  have  conducted  steady-state  numerical  analysis  in  C02  (sub-critical/super-critical)  based  loop 
with  end  heat  exchangers.  It  shows  that  if  the  average  temperature  of  the  loop  increases,  the  mass  flow  rate  and  the  heat 
transfer  rate  increases  until  the  pseudo-critical  point  and  then  decreases.  Further,  correlations  are  also  proposed  to  predict 
the  Re  and  Nu  values  in  case  of  NCL  with  end  heat  exchangers.  In  similar  study,  Yadav  et  al.  [10]  have  performed 
experimental  transient  analysis  at  sub-critical  and  supercritical  operating  conditions  (50-100  bar  and  323-353K)  by  varying 
loop  tilt  angle  and  external  mass  flow  rate.  It  was  observed  that  the  heat  transfer  rate  increases  with  increase  in  operating 
temperatures  and  external  mass  flow  rate,  whereas  it  decreases  with  the  loop  tilt  angle.  Additionally,  it  was  reported  that 
the  loop  is  quite  stable  under  the  influence  of  any  disturbances  inflected. 

Doganay  and  Turgut  [11]  studied  the  impact  of  input  power,  inclination  and  the  condenser  temperature  on  the 
thermal  performance  of  a  mini  loop  by  varying  the  particle  concentrations  of  A1203-DIW  (De-ionised  water)  nanofluid.  As 
the  concentration  of  the  nanofluids  and  inclination  increases,  the  thermal  effectiveness  of  the  system  improves.  It  was  also 
reported  that  A1203-DIW  nanofluid  has  enormous  potential  to  overhaul  the  thermal  management  problems  in  large  heat 
dissipation  devices  due  to  its  better  effectiveness  in  high-temperature  conditions.  Cheng  et  al.  [12]  have  formulated  1-D 
mathematical  model  for  loop  with  end  heat  exchangers.  The  formula  was  derived  for  temperature  distributions  along  the 
loop.  It  was  observed  that  the  modified  Grashoff  number,  heat  transfer  rate  and  Reynolds  number  increase  with  the 
increase  in  driven  temperature  difference.  Further,  Cheng  et  al.  [13]  have  theoretically  and  experimentally  studied  the  heat 
transfer  and  fluid  flow  characteristics  of  loop  with  double  pipe  heat  exchangers.  The  cooling  fluid  temperature  is  fixed 
(10°C),  and  the  heating  fluid  inlet  temperature  was  varied  from  30°C  to  60°C.  It  was  observed  that  stability  reached  for  all 
the  cases  but  the  stabilizing  time  decreases  with  an  increase  in  heating  fluid  inlet  temperature.  Upon  using  the  Nu 
correlation  for  laminar  flow,  i.e.  Nu=3.66,  significant  error  between  the  theoretical  and  experimental  results  was  observed. 
Hence  a  new  correlation  (Equation  2)  was  formulated  which  showed  much  better  agreement. 

Nu =0.2702  Re0'4713  (2) 

A  lot  of  research  highlights  analyses  in  NCL  [4-13].  However,  analysis  of  HHHC  type  is  critical  as  they 
showcase  a  twofold  instability  behaviour.  Vijayan  et  al.  [14]  have  discussed  the  stability  of  loop  with  different  heater  and 
condenser  orientations  and  experimentally  found  VHVC  and  HHHC  orientation  as  the  most  stable  and  least  stable, 
respectively.  Kudariyawar  et  al.  [15]  have  done  the  3-D  steady  and  Transient  state  characteristics  of  various  configurations 
of  NCL  viz.  VHVC,  VHHC,  HHVC  and  HHHC,  and  compared  with  experimental  data.  The  importance  of  a  low  Reynolds 
number  model  for  CFD  calculations  and  flow  initiation  transients  were  discussed.  Additionally,  the  bi-directional  and 
unidirectional  pulsations  in  the  HHHC  configuration  was  extensively  explained. 
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Many  researchers  have  discussed  the  instability  involved  in  NCL  and  techniques  to  curb  it.  Vijayan  et  al.  [1]  had 
rigorous  theoretical  and  experimental  analysis  to  prove  that  stability  increases  as  the  diameter  of  the  loop  decreases. 
Krishnani  and  Basu  [16]  performed  a  computational  study  to  see  the  influence  of  inclination  angle  on  the  stability 
behaviour  of  HHHC-type  loop.  It  was  observed  that  with  the  increase  in  power  and  sink  temperature,  the  loop  witnesses  a 
large  number  of  unstable  oscillations  and  flow  reversals.  Providing  the  loop  with  an  inclination  decreases  these  unstable 
oscillations  drastically.  The  explanation  lies  on  the  basis  that  upon  inclination,  the  buoyancy  force  reduces  by  a  factor  of 
the  cosine  of  the  angle  of  oscillation.  However,  as  inclination  increases,  the  steady-state  mass  flow  rate  of  the  system 
decreases  and  therefore,  an  ideal  inclination  of  15°  was  recommended  for  further  designs  of  NCL  systems.  Instability  is 
very  prominent  in  NCL,  especially  in  HHHC  arrangement.  Till  date  literature  reflects  various  techniques  to  overcome  this 
instability  and  uncertainty  viz.  varying  loop  diameter  [1],  presence  of  orifice  [17]  and  use  of  nanofluids  [18].  All 
mentioned  instability  restraining  techniques  have  their  own  limitations,  and  hence  the  possibility  of  Tesla  valve  (TV)  in 
loop  has  been  discussed  here  as  it  is  one  of  the  most  promising  adaptations  due  to  meagre  maintenance  requirements  and 
uni-directional  flow. 

Tesla  [2]  invented  a  valvular  conduit  to  be  used  for  unidirectional  fluid  flow.  He  stated  that  mechanical  valves 
are  not  reliable,  and  they  add  massively  to  the  cost  of  manufacture  and  maintenance.  The  basic  principle  of  the  valvular 
conduit  is  that  in  the  forward  direction,  the  fluid  flows  without  much  resistance,  but  there  is  a  significant  amount  of 
resistance  incorporated  in  the  fluid  in  the  opposite  direction  and  hence  limiting  the  flow  in  the  forward  direction  only. 
Forster  et  al.  [19]  discussed  the  design,  fabrication  and  the  testing  of  no  moving  parts  valve.  The  CFD  analysis  based 
design  was  fabricated  and  tested  positively.  Truong  and  Nguyen  [3]  had  objective  to  increase  the  value  of  diodicity  of 
TV  and  hence  optimized  the  design.  It  is  found  through  flow  simulation  that  diodicity  vary  non-linearly  with  flow  rate. 
Gamboa  et  al.  [20]  have  done  TV  shape  optimization  which  resulted  in  significant  increase  in  diodicity  values  for 
Reynolds  number  (0-2000).  The  experimental  models  were  made  using  CNC  and  verified  with  CFD  results.  The 
numerical  analysis  by  Zhang  et  al.  [21]  revealed  that  for  Reynolds  number  up  to  2000  if  the  hydraulic  diameter  is  kept 
constant,  the  TV  valve  with  a  higher  aspect  ratio  can  yield  better  diodicity.  Thompson  et  al.  [22]  have  integrated  Tesla- 
type  check  valve  to  increase  the  thermal  performance  of  a  flat-  plate  oscillating  heat  pipe  (FP-OHP).  When  it  was 
introduced,  the  desired  direction  was  promoted  with  increment  in  heat  input.  The  thermal  resistance  was  always  lesser 
in  a  TV  FP-OHP  compared  to  the  FP-OHP  without  the  valve.  Nobakht  et  al.  [23]  did  numerical  study  on  the  diodicity 
mechanism  of  different  Tesla-type  microvalves  and  found  that  the  major  pressure  loss  happens  during  the  inlet, 
irrespective  of  the  direction.  Hence  to  increase  the  diodicity,  the  flow  losses  in  the  forward  direction  has  to  be  reduced. 
The  experimental  results  of  the  microvalves  were  compared  with  numerical  simulation  and  found  good  agreement. 
Thompson  et  al.  [24]  discussed  the  transitional  and  turbulent  flow  modelling  of  TV.  The  aim  was  to  find  the  appropriate 
turbulent  model  which  agrees  the  best  with  experimental  data  and  to  optimize  the  valve  using  the  right  model  so  as  to 
increase  the  diodicity.  Further,  Thompson  et  al.  [25]  have  done  numerical  analysis  of  multi-staged  TV  and  found  that 
the  diodicity  can  increase  up  to  2  on  having  a  series  of  TVs.  This  effect  is  prominent  at  higher  Reynolds  number.  Vries 
et  al.  [26]  have  discussed  the  design  and  operation  of  TV  for  pulsating  heat  pipes.  It  was  designed  on  the  basis  of  2-D 
and  3-D  numerical  analysis  for  enhancing  the  flow  in  a  pulsating  heat  pipe.  The  aim  was  to  maximize  the  diodicity  and 
retain  low-pressure  drop  in  the  promoted  direction.  In  the  numerical  analysis,  the  diodicity  found  was  higher  for  higher 
Reynolds  numbers.  Even  Rajat  [27]  has  mentioned  the  benefit  of  adding  TV  in  NCL  as  numerical  result  gave 
unidirectional  flow. 
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Figure  4:  Fabricated  TV.  Figure  5:  NCL  with  TV. 


Figure  6:  Experimental  set  up. 


To  test  the  stability  of  NCL,  power  variation  made  was  90W,  110W,  130W,  150W  and  170W.  Additionally,  the 
loop  was  inclined  at  angles  viz.  2.5°,  5°  and  10°  to  understand  the  working  of  the  TV  at  these  angles  and  to  compare  them 
without  the  TV. 


5.  RESULTS  AND  DISCUSSIONS 


Here,  the  main  parameters  for  data  reduction  are  temperatures  at  different  points  and  pressure  within  the  loop.  The  stability 
of  TV-based  loop  was  tested  in  two  ways,  firstly  by  measuring  the  differential  pressure  between  two  legs  and  secondly  by 
analyzing  the  variation  and  difference  in  temperature  at  the  lower  side  of  each  leg  (lower  right  -  lower  left).  Pressure  curve 
shows  stability  when  the  oscillations  are  reduced  to  minimum.  A  constant  non-zero  value  of  AT  signifies  the  formation  of  a 
loop,  i.e.  counterclockwise  (positive)  and  clockwise  (negative)  when  viewing  from  the  front.  The  tests  covered  heater 
power  ranging  from  90W  to  170W  with  a  step  of  20W. 

The  temperature  variations  give  insights  into  the  behaviour  of  the  loop  over  time  at  different  points  in  the  loop. 
For  a  heater  power  of  150  W,  it  can  be  seen  from  Figure  7  that  the  temperature  profile  of  the  loop  with  TV  is  smoother 
than  the  one  without  TV  (Figure  8).  The  time  taken  to  reach  stability  with  TV  is  also  less  with  respect  to  the  loop  without 
TV. 

Further,  in  the  case  of  TV-based  loop,  the  AT  graph  is  positive,  indicating  fluid  movement  in  forward  direction  of 
the  Tesla.  As  mentioned  earlier,  AT  graph  is  smoother  with  TV  for  different  heater  input  power.  At  90  W  power,  the  loop 
with  TV  has  reached  more  stable  state  faster  than  that  of  the  loop  without  TV,  which  shows  the  stabilizing  effect  TV  has  on 
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NCL.  Also,  the  high  and  steep  peaks  of  the  AT  and  AP  of  the  TV  during  the  initial  stages  of  heating  reveals  that  the  TV  is 
more  responsive  to  the  heating  effects  of  the  fluid  as  shown  in  Figure  9.  A  similar  trend  can  be  observed  in  the  AP  values 
for  loop  without  TV  which  is  close  to  zero  as  shown  in  Figure  10.  It  means  the  dynamic  pressure  difference  in  the  loop  is 
negligible,  and  hence  the  flow  is  minimal  or  erratic.  In  the  loop  with  TV,  the  pressure  difference  has  a  non-zero  value, 
signifying  the  presence  of  noticeable  amount  of  flow  and  that  too  in  forward  direction.  Similarly  at  170  W  also  the  instability  in 
terms  of  variations  in  temperature  (Figure  11)  and  pressure  (Figure  12)  for  loop  without  TV  which  takes  more  time  to  attain  a 
constant  value  compared  to  the  case  of  TV. 

Next,  AT  and  AP  plots  are  made  considering  all  power  trials  (90  W-170W).  They  show  a  similar  trend,  having  an 
initial  peak  and  then  after  a  couple  of  oscillations  reaches  more  stable  temperature  difference.  The  general  observation  made 
here  is  the  value  of  temperature  difference  which  is  more  when  the  wattage  is  increased.  In  the  beginning,  the  peak  in  AT  and  AP 
occurs  due  to  initial  heating  in  the  horizontal  heater  (which  can  take  any  direction).  Further,  in  the  loop  without  TV  it  can  be 
either  of  the  direction  and  in  the  loop  with  TV,  after  initial  heating,  the  flow  is  biased  towards  forward  direction  always  and  in 
smoother  way  as  shown  in  Figure  13.  The  peak  observed  in  AP  is  also  due  to  the  sudden  rising  of  the  hot  fluid  from  the  lower 
portion  to  the  upper  portion,  but  as  seen  from  Figure  14,  TV-based  loop  performs  better. 


Figure  7:  Temperature  Profile  with  TV.  Figure  8:  Temperature  Profile  without  TV. 


Figure  9:  AT  comparison  for  90W.  Figure  10:  AP  comparison  for  90W. 
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Figure  11:  AT  Comparison  for  170  W. 


Figure  12:  AP  Comparison  for  170  W. 


(a)  Loop  with  Tesla  valve.  (b)  Loop  without  Tesla  valve. 


Figure  13:  AT  vs  Time 


(a)  Loop  with  Tesla  valve.  (b)  Loop  without  Tesla  valve. 

Figure  14:  AP  vs  Time. 
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Figure  15:  AT  at  Different  Inclinations  with  TV.  Figure  16:  AP  at  Different  Inclinations  without  TV. 

Even,  AP  and  AT  curves  are  also  analyzed  to  understand  the  flow  behaviour  when  the  loop  was  inclined.  In  the 
loop  with  Tesla  valve,  the  forward  direction  of  the  fluid  is  indicated  by  the  plot  with  positive  AT.  AP  is  used  as  an 
indication  of  the  amount  of  flow  within  the  loop.  The  effect  of  TV  is  tested  at  0°,  2,5°  and  5°  inclinations  (loop  was 
inclined  so  as  to  oppose  the  effect  of  TV  and  prefered  the  flow  in  the  reverse  direction)  with  150W  heater  power. 
According  to  the  AT  graph  (Figure  15),  the  flow  with  TV  at  0°  inclination  took  forward  direction  (AT  is  positive),  and  the 
flow  in  the  loop  at  other  angles  took  reverse  (AT  is  negative).  This  shows  that  TV  is  not  able  to  overcome  the  flow  bias 
offered  due  to  the  inclination  of  the  loop  and  works  better  at  horizontal  position.  From  AP  graph  (Figure  16),  the  highest 
flow  resistance  at  loop  angle  2.5°  is  observed,  which  indicates  the  effect  of  TV.  Additionally,  it  can’t  resist  the  flow  at 
angle  more  than  2.5°. 

6.  CONCLUSIONS 

Effect  of  TV  in  NCL  was  studied  by  using  two  separate  loops  subjected  to  the  same  test  conditions.  Some  of  the  critical 
observations  made  out  of  the  experimental  investigation  are  summarised  below. 

•  Directionality:  For  the  entire  range  of  heater  power,  the  loop  with  TV  showed  uni-directional  behaviour  and 
favoured  the  flow  in  forward  direction. 

•  Stability  Time:  From  the  graphs,  it  was  observed  that  the  time  taken  for  the  loop  to  reach  stability  was  lower  in 
the  TV-based  loop  compared  to  the  without  TV. 

•  Effect  of  Power:  The  importance  of  a  TV  can  be  better  observed  on  varying  power  input.  The  AT  curves  for  loop 
with  TV  had  similar  smooth  profile,  whereas  in  the  loop  without  TV  there  is  no  observable  pattern,  indicating 
lesser  control  on  the  flow  inside  the  loop. 

•  Effect  of  Inclination:  In  inclined  position,  the  loop  without  TV  readily  took  the  direction  favoured  by  the 
inclination.  However,  the  loop  with  TV  offered  considerable  resistance  to  the  flow  until  2.5°  angle  of  inclination. 
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